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Abstract 

We propose a general geometric growth model for pseudofractal scale-free web, 
which is controlled by two tunable parameters. We derive exactly the main charac- 
teristics of the networks: degree distribution, second moment of degree distribution, 
degree correlations, distribution of clustering coefficient, as well as the diameter, 
which are partially determined by the parameters. Analytical results show that the 
resulting networks are disassortative and follow power-law degree distributions, with 
a more general degree exponent tuned from 2 to 1 + ; the clustering coefficient of 
each individual node is inversely proportional to its degree and the average cluster- 
ing coefficient of all nodes approaches to a large nonzero value in the infinite network 
order; the diameter grows logarithmically with the number of network nodes. All 
these reveal that the networks described by our model have small-world effect and 
scale-free topology. 
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1 Introduction 



Since the pioneering papers by Watts and Strogatz on small- world networks [1] 
and Barabasi and Albert on scale- free networks [2], complex networks, which 
describe many systems in nature and society, have become an area of tremen- 
dous recent interest [3,4,5,6,7]. In the last few years, modeling real-life systems 
has attracted an exceptional amount of attention within the physics commu- 
nity. While a lot of models have been proposed, most of them are stochastic 
[3,4,5,6,7]. However, because of their advantages, deterministic networks have 
also received much attention [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31]. 
First, the method of generating deterministic networks makes it easier to gain 
a visual understanding of how networks are shaped, and how do different nodes 
relate to each other [8]; moreover, deterministic networks allow to compute an- 
alytically their properties: degree distribution, clustering coefficient, average 
path length, diameter, betweenness, modularity and adjacency matrix whose 

eigenvalue spectrum characterizes the topology [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,2 

The first model for deterministic scale-free networks was proposed by Barabasi 
et al. in Ref. [8] and was intensively studied in Ref. [9]. Another elegant model, 
called pseudofractal scale-free web (PSW) [10], was introduced by Dorogovt- 
sev, Goltsev, and Mendes, and was extended by Cornelias et al. in Ref. [11]. 
Based on a similar idea of PSW, Jung et al. presented a class of recursive 
trees [12]. Additionally, in order to discuss modularity, Ravasz et al. proposed 
a hierarchical network model [13,14], the exact scaling properties and extensive 
study of which were reported in Refs. [15] and [16], respectively. Recently, In re- 
lation to the problem of Apollonian space-filing packing, Andrade et al. intro- 
duced Apollonian networks [17] which were also proposed by Doye and Massen 
in Ref. [18] and have been intensively investigated [22,23,32,33,34,35,36]. In 
addition to the above models, deterministic networks can be created by vari- 
ous techniques: modification of some regular graphs [24], addition and product 
of graphs [25], edge iterations [26,27] and other mathematical methods as in 
Refs. [28,29,30,31]. 

As mentioned by Barabasi et al., it would be of major theoretical interest 
to construct deterministic models that lead to scale-free networks [8]. Here 
we do an extensive study on pseudofractal scale-free web [10]. The PSW can 
be considered as a process of edge multiplication. In fact, a clique (edge is 
a special case of it) can also reproduce new cliques and the number of the 
new reproduction may be different at a time. Motivated by this, in a simple 
recursive way we propose a general model for PSW by including two param- 
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Fig. 1. Scheme of the growth of the deterministic network for the case of m = 2 and 
q = 2. Only the first three steps are shown. 

eters, with PSW as a particular case of the present model. The deterministic 
construction of our model enables one to obtain the analytic solutions for its 
structure properties. By adjusting the parameters, we can obtain a variety of 
scale-free networks. 



2 The general geometric growth model 

Before introducing our model we give the following definitions on a graph 
(network). The term size refers to the number of edges in a graph. The num- 
ber of nodes in a graph is called its order. When two nodes of a graph are 
connected by an edge, these nodes are said to be adjacent, and the edge is said 
to join them. A complete graph is a graph in which all nodes are adjacent to 
one another. Thus, in a complete graph, every possible edge is present. The 
complete graph with q nodes is denoted as K q (also referred in the literature 
as q- clique). Two graphs are isomorphic when the nodes of one can be rela- 
beled to match the nodes of the other in a way that preserves adjacency. So 
all g-cliques are isomorphic to one another. 

The network is constructed in a recursive way. We denote the network after 
t steps by R(q,t), q > 2,t > (see Fig. 1). Then the network at step t is 
constructed as follows: For t = 0, R(q, 0) is a complete graph K q+1 (or (q+ 1)- 
clique) consist of q + 1 g-cliques), and R(q, 0) has q + 1 nodes and q(q + l)/2 
edges. For t > 1, R(q,t) is obtained from R(q,t — 1) by adding m new nodes 
for each of its existing subgraphs isomorphic to a g-clique, and each new node 
is connected to all the nodes of this subgraph. In the special case m = 1 and 
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Fig. 2. Growth process of pseudofractal scale-free web (i.e. the case of m = 1 and 
g = 2), exhibiting the first three steps. 



Fig. 3. Aggregation process from R(2, 1) to R(2, 2), which is obtained by adjoining 
of three copies of R(2, 1) at the hubs. 

q — 2, it is reduced to the pseudofractal scale-free web described in Ref. [10]. 
In the limiting case of m — 1, we obtain the same networks as in Ref. [11]. 
However, our family is richer as m can take any natural value. 

There is an interpretation called 'aggregation' [37] for our model. As an ex- 
ample, here we only explain them for the case of m — 1 and q = 2. Figure 
2 illustrates the growing process for this particular case, which may be ac- 
counted for as an 'aggregation' process described in detail as follows. First, 
three of the initial triangle (t = 0) are assembled to form a new unit (t = 1). 
Then we assemble three of these units at the hubs (the nodes with highest 
degree) in precise analogy with the step leading from t — to t — 1 to form a 
new cell (t = 2) (see Fig. 3). This process can be iterated an arbitrary number 
of times. Moreover, an alternative explanation of our model which is often 
useful is that of 'miniaturization' (see Ref. [37]). 



3 Structural properties 

Below we will find that the tunable parameters m and q control some relevant 
characteristics of the network R(q,t). Because q = 2 is a particular case, for 
conveniences, we treat q = 2 and q > 3 separately. 
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3.1 Limiting case of q = 2 



Order and size. In the case of q = 2, we denote R(q, t) by Q{t). Let us consider 
the total number of nodes N t and total number of edges E t in Q(t). Denote 
n v (t) as the number of nodes created at step t. Note that the addition of each 
new node leads to two new edges. By construction, for t > 1, we have 

n v (t) = mE t -i (1) 



and 

E t = E t _ 1 + 2n v (t). (2) 

Considering the initial condition n v (0) = 3 and E = 3, it follows that 

n v (t) = 3m(2m+ l)*' 1 and E t = 3 (2m + 1)*. (3) 



Then the number of nodes increases with time exponentially and the total 
number of nodes present at step t is 



*r ,+ s 3 2m + l* + 3 

N t =J2Mti) = - ^ • (4) 

*i=0 z 

Thus for large t, The average degree k t = ^ is approximately 4. 

Degree distribution. Let fcj(i) be the degree of node i at step t. Then by con- 
struction, it is not difficult to find following relation: 

ki(t) = (m + l)ki(t-l), (5) 



which expresses a preference attachment [2] . If node % is added to the network 
at step ti, ki(ti) = 2 and hence 

fc i (i) = 2(m + l) t - ti . (6) 
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Therefore, the degree spectrum of the network is discrete. It follows that the 
degree distribution is given by 



P(k) = 



f Mo) 

N t 
n v {U) 

N t 





(2m+l)* + l 
2m(2m + I)*'- 1 

(2m +!)* + ! 



for h = 



for ti > 1 



otherwise 



(7) 



and that the cumulative degree distribution [6,10] is 
MP) (2m + l)' l + l 



(2m +!)* + ! 



(8) 



Substituting for ti in this expression using ti = t — ^(m+i) gives 



-Pcum(^) 



_ ln(2m+l) 

(2m+l)*(fc/2) m™+d + 1 



(2m+ 1)* + 1 

ln(2m+l) 



fc\ ln(m+l) 



for large t. 



(9) 



So the degree distribution follows the power law with the exponent 7 = 1 + 
^nim+rj ■ ^ or ^ e particular case of m = 1, Eq. (9) recovers the result previously 
obtained in Ref. [10]. 

Second moment of degree distribution. Let us calculate the second moment of 
degree distribution (A; 2 ). It is defined by 



(k 2 ) = ^i2n v ( tl )[k(t u t)}\ (10) 
iV * ti=0 

where k(ti, t) is the degree of a node at step t, which was generated at step tj. 
This quality expresses the average of degree square over all nodes in the net- 
work. It has large impact on the dynamics of spreading [38,39] and the onset of 
percolation transitions [40,41] taking place in networks, when {k 2 } is diverging, 
the networks allow the onset of large epidemics whatever the spreading rate 
of the infection [38,39], at the same time the networks are extremely robust to 
random damages, in other words, the percolation transition is absent [40,41]. 

substituting Eqs. (3), (4) and (6) into Eq. (10), we derive 
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(k 2 



1 + 1 

1 ^ (2m+l)' 

8 (m + l) 2f+1 
m(2m + 1) 




1 

m 



In this way, second moment of degree distribution (k 2 ) has been calculated 
explicitly, and result shows that it diverges as an exponential law. So the 
networks are resilient to random damage and are simultaneously sensitive to 
the spread of infections. 

Degree correlations. As the field has progressed, degree correlation [42,43,44,45,46] 
has been the subject of particular interest, because it can give rise to some 
interesting network structure effects. An interesting quantity related to degree 
correlations is the average degree of the nearest neighbors for nodes with de- 
gree k, denoted as k nn (k) [43,44]. When k mi (k) increases with k, it means that 
nodes have a tendency to connect to nodes with a similar or larger degree. In 
this case the network is defined as assortative [45,46]. In contrast, if k nn (k) is 
decreasing with k, which implies that nodes of large degree are likely to have 
near neighbors with small degree, then the network is said to be disassortative. 
If correlations are absent, k nn (k) = const. 

We can exactly calculate k nn for the networks using Eq. (6) to work out how 
many links are made at a particular step to nodes with a particular degree. 
Except for three initial nodes generated at step 0, no nodes born in the same 
step, which have the same degree, will be linked to each other. All links to 
nodes with larger degree are made at the creation step, and then links to 
nodes with smaller degree are made at each subsequent steps. This results in 
the expression 



*;=* \ 

m-n v (t i )k(t i ,t' i -l)k(t' i ,t)\ (12) 

for A; = 2 (m + 1)*~*\ Here the first sum on the right-hand side accounts 
for the links made to nodes with larger degree (i.e. t\ < U) when the node 
was generated at U. The second sum describes the links made to the current 
smallest degree nodes at each step t\ > ti. 

Substituting Eqs. (3) and (6) into Eq. (12), after some algebraic manipulations, 
Eq. (12) is simplified to 
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2(2m + 1) 



m 



(m + If 



2(m + l) 2m . 

+ —(t-U). 



m 



m + 1 



2m + 1 

Thus after the initial step k nn grows linearly with time. 
Writing Eq. (13) in terms of k, it is straightforward to obtain 



(13) 



2(2m + 1) 



m 



(m + 1) 



2m + 1 
2(m+l) 




(m+ir 

2m+l 



i(m + l) 



+ 



ln(|) 



m 



m + 1 ln(m + 1) 



(14) 



Therefore, k nn (k) is approximately a power law function of k with negative 
exponent, which shows that the networks are disassortative. Note that k nn (k) 
of the Internet exhibit a similar power-law dependence on the degree k nn (k) ~ 
AT", with to = 0.5 [43]. 

Clustering coefficient. The clustering coefficient defines a measure of the level 
of cohesiveness around any given node. By definition, the clustering coeffi- 
cient [1] Cj of node i is the ratio between the number of edges e, that actu- 
ally exist among the fcj neighbors of node i and its maximum possible value, 
ki(ki — l)/2, i.e., C% = 2ei/ki(ki — l). The clustering coefficient of the whole net- 
work is the average of all individual C-s. Next we will compute the clustering 
coefficient of every node and their average value. 

Obviously, when a new node % joins the network, its degree kj and is 2 and 
1, respectively. Each subsequent addition of a link to that node increases both 
hi and by one. Thus, equals to fcj — 1 for all nodes at all steps. So one 
can see that, there is a one-to-one correspondence between the degree of a 
node and its clustering. For a node with degree k, the exact expression for 
its clustering coefficient is 2/k. Therefore, the clustering coefficient spectrum 
of nodes is discrete. Using this discreteness, it is convenient to work with the 
cumulative distribution of clustering coefficient [10] as 



W cnm (C) 



-E 



2n v (tj) 
k(tj,t) 



lnfgm+l) 

C ln ( m+1 ) = C 7-1 . 



(15) 



It is worth noting that for the special case of m = 1, this result has been 
obtained previously [10]. 

The clustering coefficient of the whole network at arbitrary step t can be easily 
computed, 
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1 « 2 n v (r) _ 2m + 2 (m + l) f (2m + 1)* + g±f 
* N t f^ Q k{r,t) 2m + 3 (m + l)*[(2m+ 1)* + 1] ' 1 J 

In the infinite network size limit (i — > oo), C = {2m + 2) /{2m + 3). Thus 
the clustering C is high and increases with m. Moreover, similarly to the 
degree exponent 7, C is tunable by choosing the right value of parameter m: 
in particular, C ranges from 4/5 (in the special case of m — 1 [10]) to limit of 
1 when m becomes very large. 

Diameter. The diameter of a network is defined as the maximum of the shortest 
distances between all pairs of nodes, which characterizes the longest commu- 
nication delay in the network. Small diameter is consistent with the concept 
of small-world and it is easy to compute for our networks. Below we give the 
precise analytical computation of diameter of Q(t) denoted by Diam{Q{t)). 

It is easy to see that at step t — (resp. t — 1), the diameter is equal to 1 
(resp. 2). At each step t > 2, one can easily see that the diameter always lies 
between a pair of nodes that have just been created at this step. In order to 
simplify the analysis, we first note that it is unnecessary to look at all the 
nodes in the networks in order to find the diameter. In other words, some 
nodes added at a given step can be ignored, because they do not increase the 
diameter from the previous step. These nodes are those that connect to edges 
that already existed before step t — 1. Indeed, for these nodes we know that a 
similar construction has been done in previous steps, so we can ignore them for 
the computation of the diameter. Let us call "outer" nodes the nodes which 
are connected to a edge that did not exist at previous steps. Clearly, at each 
step, the diameter depends on the distances between outer nodes. 

At any step t > 2, we note that an outer node cannot be connected with two 
or more nodes that were created during the same step t > t — 1. Indeed, we 
know that from step 2, no outer node is connected to two nodes of the initial 
triangle &(0). Thus, for any step t > 2, any outer node is connected with 
nodes that appeared at pairwise different steps. Now consider two outer nodes 
created at step t > 2, say v t and w t . Then v t is connected to two nodes, and 
one of them must have been created before or during step t — 2. We repeat this 
argument, and we end up with two cases: (1) t = 2m is even. Then, if we make 
m "jumps", from Vt we reach the initial triangle G(0), in which we can reach 
any w t by using an edge of C?(0) and making m jumps to w t in a similar way. 
Thus Diam(G(2m)) < 2m + 1. (2) t = 2m + 1 is odd. In this case we can stop 
after m jumps at (?(1), for which we know that the diameter is 2, and make m 
jumps in a similar way to reach w t . Thus Diam(Q(2m + 1)) < 2{m + 1). It is 
easily seen that the bound can be reached by pairs of outer nodes created at 
step t. More precisely, those two nodes v t and w t share the property that they 
are connected to two nodes that appeared respectively at steps t — 1, t — 2. 
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Hence, formally, Diam(Q(t)) — t + 1 for any t > 0. Note that iV* ~ (2m + 
1)*, thus the diameter is small and scales logarithmically with the number of 
network nodes. 



3.2 Case of q> 3 



In these cases, the analysis is a little difficult than those of the last subsection. 
An alternative approach has to be adopted, although it may also holds true for 
the first case in some situations. The method of the last subsection is relatively 
easy to generalize to these cases, and below we will address it, focusing on 
order, size, degree distribution, clustering coefficient and diameter. 

Order and size. Let n v (t), n e (t) be the number of nodes and edges created 
at step t, respectively. Denote K qit as the total number of g-cliques in the 
whole network at step t. Note that the addition of each new node leads to 
q new g-cliques and q new edges. By construction, we have n e (t) = qn v (t), 
n v (t) = mK qt _x and K qt = K qt _ x + qn v (t). Thus one can easily obtain 
K q , t = (mq + ^Kg^ = (g + l)(mg + l)< (t > 0), n v (t) = m(q + l)(mq + 
(t > 0) and n e (t) = mq(q+ l)(mq + (t > 0). From above results, we can 
easily compute the order and size of the networks. The total number of nodes 
N t and edges E t present at step t is 



N t = E n v (U) 



(q + l)[(mq + iy + q-l\ 



(17) 



and 



Et = E n e (U) = (q + l)(mq + 1)* + (g + 2) , (18) 

respectively. For infinite t, the average degree k t = ^ is approximately 2q. 

Degree distribution. When a new node % is added to the graph at step ti, it 
has degree q and forms q new g-cliques. Let n q (i,t) be the total number of 
g-cliques at step t that will created new nodes connected to the node i at 
step t + 1. So at step U, n q (i,ti) = q. By construction, we can see that in the 
subsequent steps each new neighbor of i generates q — 1 new g-cliques with i 
as one node of them. Let ki(t) be the degree of % at step t. It is not difficult 
to find following relations for t > ti + 1: 

Afci(t) = hit) -ki(t-l)= mn q (i, t - 1) (19) 
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and 

n q (i, t) = n q (i, t - 1) + (q - l)A^(t). (20) 

From the above two equations, we can derive n q (i, t) = [m(q — l) + l]n q (i, t — 1). 
Considering n q (i, U) = q, we obtain n q (i, t) = q[m(q — 1) + l]* - *' and Aki(t) = 
mq[m(q — 1) + Then the degree ki(t) of node % at time t is 



t 

k i (t)=k i (t i )+ Ak i( f h) 

t h =U+l 

_ q[m(q - 1) + l} 1 ' 1 ' + q 2 - 2q 



(21) 



Since the degree of each node has been obtained explicitly as in Eq. (21), 
we can get the degree distribution via its cumulative distribution [6,10], i.e. 
Pcum(k) = Y,k'>kN(k',t)/N t ~ A; 1-7 , where N(k',t) denotes the number of 
nodes with degree k' . The analytic computation details are given as follows. 
For a degree k 

_ q[m(q - 1) + + q 2 - 2q 
k - — , (22) 



there are n v (j) = m(q + l)(qm + l)- 7 ' -1 nodes with this exact degree, all of 
which were born at step j. All nodes with birth time at j or earlier have this 
and a higher degree. So we have 

£ *(*-, t) = ± n„( S ) = fa+DK"* + 1)^-1] . (23) 



As the total number of nodes at step t is given in Eq. (17), we have 



q[m(q — 1) + 1J 3 + q — 2q 



q-l 



{q+l)[(rnq+iy+q-l\ 

q 



(24) 



Therefore, for large t we obtain 

[m(q - 1) + 1_7 = (mq + 1) J_ * (25) 



and 



1SS 1 + Ere^Tif (26) 
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Fig. 4. The dependence relation of C on q and m. 

For the special case m — 1, Eq. (26) recovers the results previously reported 
in Ref. [11]. 

Clustering coefficient. The analytical expression for clustering coefficient C(k) 
of the individual node with degree k can be derived exactly. When a node is 
created it is connected to all the nodes of a g-clique whose nodes are completely 
interconnected. Its degree and clustering coefficient are q and 1, respectively. 
In the following steps, if its degree increases one by a newly created node 
connecting to it, then there must be q — 1 existing neighbors of it attaching 
to the new node at the same time. Thus for a node of degree k, we have 



C(k) = 



+ (q-l)(k-q) _ 2(g-l)(fc- 
W-V fc(jfe-l) 



(27) 



which depends on degree k and q. For k ^> q, the C(fc) is inversely proportional 
to node degree. The scaling C(k) ~ has been found for some network 
models [4,5,6,7], and has also observed in several real-life networks [14]. 

Using Eq. (27), we can obtain the clustering C t of the networks at step t: 
Ct_ A^o D r .(D r -l) ' (28) 



where the sum is the total of clustering coefficient for all nodes and D r = 
q [ m (g-i)+i] t --W-2q shown by Eq> (2i) i s the degree of the nodes created at 
step r. 

It can be easily proved that for arbitrary fixed m, C t increases with q, and that 
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for arbitrary fixed g, Ct increases with m. In the infinite network order limit 
(N t — > oo), Eq. (28) converges to a nonzero value C. When q = 2, for m — 1, 2, 
3 and 4, C equal to 0.8000, 0.8571 0.8889 and 0.9091, respectively. When m = 
2, for q = 2, 3, 4 and 5, C are 0.8571, 0.9100, 0.9348 and 0.9490, respectively. 
Therefore, the clustering coefficient of our networks is very high. Moreover, 
similarly to the degree exponent 7, clustering coefficient C is determined by 
q and m. Figure 4 shows the dependence of C on q and m. 

Diameter. In what follows, the notations \x\ and |_^J express the integers ob- 
tained by rounding x to the nearest integers towards infinity and minus infinity, 
respectively. Now we compute the diameter of R(q,t), denoted Diam(R(q,t)) 
for q > 3 (q = 2 is a particular case that is treated separately in the last 
subsection) : 

Step 0. The diameter is I. 

Steps 1 to [^tt-] — 1. In this case, the diameter is 2, since any new node is by 
construction connected to a g-clique forming a (q + l)-clique, and since any 
(q + l)-clique during those steps contains at least \ q ^r\ (q even) or [ 2 |^-]+l (q 
odd) nodes from the initial (q + l)-clique R(q, 0) obtained after step 0. Hence, 
any two newly added nodes u and v will be connected respectively to sets S u 
and S v , with S u C V(R(q,0)) and S v C V(R(q,0)), where V(R(q,0)) is the 
node set of R(q,0); however, since \S U \ > f 2 ^-] (g even) and \S V \ > [^y^+l 
(q odd), where \S\ denotes the number of elements in set S, we conclude that 
S U D S v ^ 0, and thus the diameter is 2. 

Steps r^Tr-] to q. In any of these steps, some newly added nodes might not 
share a neighbor in the original (q + l)-clique R(q,0); however, any newly 
added node is connected to at least one node of the initial (q + l)-clique 
R(q,0). Thus, the diameter is equal to 3. 

Further steps. Similar to the case of q — 2, we call "outer" nodes the nodes 
which are connected to a g-clique that did not exist at previous steps. Clearly, 
at each step, the diameter depends on the distances between outer nodes. Now, 
at any step t > q + 1, an outer node cannot be connected with two or more 
nodes that were created during the same step < t' < t — 1. Moreover, by 
construction no two nodes that were created during a given step are neighbors, 
thus they cannot be part of the same g-clique. Therefore, for any step t > g+1, 
some outer nodes are connected with nodes that appeared at pairwise different 
steps. Thus, if v t denotes an outer node that was created at step t, then v t is 
connected to nodes i>jS, 1 < i < t — 1, where all the is are pairwise distinct. 
We conclude that Vt is necessarily connected to a node that was created at a 
step to < t — q. If we repeat this argument, then we obtain an upper bound 
on the distance from v t to the initial (g + l)-clique R(q, 0). Let t = aq + p, 
where 1 < p < q. Then, we see that v t is at distance at most cc + 1 from a node 
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in R(q,t). Hence any two nodes Vt and Wt in R(q,0) lie at distance at most 
2(a + 1) + 1 ; however, depending on p, this distance can be reduced by 1, since 
when p < — 1, we know that two nodes created at step p share at least a 
neighbor in R(q,0). Thus, when 1 <p < [2±1] -1, Diam(R(q,t)) < 2(a + l), 
while when |"2±I] < p < q, Diam(R(q,t)) < 2(a + 1) + 1. One can see that 
these bounds can be reached by pairs of outer nodes created at step t. More 
precisely, those two nodes v t and Wt share the property that they are connected 
to q nodes that appeared respectively at steps t — l,t — 2, . . .t — q. 

Based on the above arguments, one can easily see that for t > q, the diameter 
increases by 2 every q steps. More precisely, we have the following result, for 
any q > 3 and t > 1 (when t — 0, the diameter is clearly equal to 1): 

Diam(R(q,t)) = 2(L^J + 1) + /(?,*), 

where f(q,t) — if t — [^-\q < f 2 ^"! — 1> an d 1 otherwise. When t gets 
large, Diam(R(q,t)) ~ ^, while JV t ~ (mq + 1)*, thus the diameter grows 
logarithmically with the number of nodes. 

It is easy to see that these cases of q > 3 have very similar topological prop- 
erties to the case q — 2. Additionally, for the cases of q > 3, the networks will 
again be disassortative with respect to degree because of the lack of links be- 
tween nodes with the same degree; the second moment of degree distribution 
(k 2 ) will also diverge, which is due to the fat tail of the degree distribution. 



4 Conclusion and discussion 

To sum up, we have proposed and investigated a deterministic network model, 
which is constructed in a recursive fashion. Our model is actually a tunable 
generalization of the growing deterministic scale-free networks introduced in 
Ref. [10]. Aside from their deterministic structures, the statistical properties 
of the resulting networks are equivalent with the random models that are 
commonly used to generate scale-free networks [4,5,6,7]. We have obtained the 
exact results for degree distribution and clustering coefficient, as well as the 
diameter, which agree well with large amount of real observations [4,5,6,7]. The 
degree exponent can be adjusted, the clustering coefficient is very large, and 
the diameter is small. Therefore, out model may perform well in mimicking a 
variety of scale-free networks in real-life world. Moreover, our networks consist 
of cliques, which has been observed in variety of the real-world networks, such 
as movie actor collaboration networks, scientific collaboration networks and 
networks of company directors [4,5,6,7]. 
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